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We propose a model for an anisotropic dark energy star where we assume that the radial pressure 
exerted on the system due to the presence of dark energy is proportional to the isotropic perfect 
fluid matter density. We discuss various physical features of our model and show that the model 
satisfies all the regularity conditions and is stable as well as singularity-free. 
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I. INTRODUCTION 

Current cosmological observations of the accelerated 
expansion of the universe strongly suggest that about 
96% of the total energy content of the universe is exotic 
in nature out of which 73% is believed to be gravitation- 
ally repulsive in nature popularly called dark energy and 
the remaining 23% is attractive in nature and exists in 
the form of dark matter [l|, ■ Consequently, cosmolog- 
ical models based on dark energy either in the form of 
a cosmological constant or in some other exotic forms of 
matter have got tremendous attention in the recent past. 
From the astrophysical perspective, if it is fundamentally 
impossible to get any observational evidence for the exis- 
tence of an event horizon in our universe Q (though our 
current understanding of the general theory of relativ- 
ity strongly favour the existence of strong gravitational 
regions induced by compact objects covered under the 
event horizon), one is tempted to look for alternative 
models which may serve as alternatives to black holes. 
A dark energy star is, in particular, interesting in this 
scenario [4]. 

Any interior solution to the vacuum Schwarzschild ex- 
terior comprising a fluid distribution governed by an 
equation of state(EOS) of the form p — —\p, may be 
considered as a dark energy star !5]. In the past, vari- 
ous model specific dark energy stars have been proposed 
(see for example, Ref. [J-Q and Ref. [9] for a recent 
review). In the present work, we propose a model for an 
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anisotropic dark energy star where we assume that the ra- 
dial pressure exerted on the system due to the presence of 
dark energy is proportional to the isotropic perfect fluid 
matter density. The stellar configuration comprises two 
fluids - an ordinary baryonic perfect fluid together with 
an yet unknown form of matter (dark energy) which is 
repulsive in nature. We also assume that the two fluids 
are non interacting amongst each other. 

To describe the energy-momentum tensor for such a 
hybrid model, we assume that our resulting composition 
is anisotropic in nature, i.e. pr ^ pt, where Pr and pt cor- 
respond to radial and tangential pressure, respectively. 
Ever since the pioneering works of Bowers and Liang 
(lOj . anisotropic relativistic stellar models have played 
an inportant role in the description of compact stellar 
objects (see [ll| for a recent review). At the microscopic 
level, a variety of reasons such as the existence of type ZA 
superfluid, phase transition, pion condensation, rotation, 
magnetic fleld, mixture of two fluids, bosonic composition 
etc., may give rise to anisotropic pressures inside a stel- 
lar object. Recent observations on highly compact astro- 
physical objects like X ray pulsar Her X-1, X ray buster 
4U 1820-30, millisecond pulsar SAX J 1808.4 - 3658, X 
ray sources 4U 1728 - 34, etc., also strongly favour an 
anisotropic matter distribution since the density inside 
such an ultra-compact object is expected to be beyond 
nuclear matter density. In our model, we assume that 
the ansisotropy is generated due to two kinds of fluid 
distributions. In a repent paper [T2j . an exterior solution 
corresponding to a two fluid stellar model composed of 
non-interacting phantom scalar field describing the dark 
energy and ordinary matter has been reported which re- 
duces to Schwarzschild solution in the absence of the dark 
energy. In our paper, we match the interior solution to 
the Schwarzschild exterior solution at the boundary. 

Our paper is organized as follows: In Section II we 
have provided the basic equations in connection to the 
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proposed model for dark energy star. Sections. Ill - 
VIII are dealt, respectively, with the boundary condi- 
tions, TOV equation, energy conditions, stability, mass- 
radius relation and junction conditions for the solutions 
under consideration. Some concluding remarks are made 
in the Section IX. 



II. BASIC EQUATIONS AND THEIR 
SOLUTIONS 



To describe the space-time of the dark en ergy stellar 
configuration, we take the Krori and Barua [13| metric 
(henceforth KB) given by 

ds^ = -^e-^'^'-^dt^ + e^'-^'Ur^ + r^d9^ + sin'^Bd(t>'^), (1) 

with A(r) = Ar^ and v{r) = Br^ + C where A, B and 
C are arbitrary constants to be determined on physical 
grounds. The energy-momentum of the two fluids are 
such that 
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where p and p correspond to the energy density and pres- 
sure of the baryonic matter, respectively, and pde, Pder 
and pdet are the 'dark' energy density, radial pressure and 
tangential pressure, respectively. The left hand sides of 
equations ©-(HI) are the effective energy- density and two 
pressures, respectively, of the composition. 

The Einstein's field equations for the metric ([T]) are 
then obtained as (we assume G = c — 1 under ge- 
ometrized relativistic units) 
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To solve the above set of equations, we assume that 
the dark energy radial pressure is proportional to the 
dark energy density, i.e.. 



Pde 



-pde 



(8) 



and the dark energy density is proportional to the matter 
density, i.e.. 



Pde 



ap, 



(9) 



where a > is a proportionality constant. In connec- 
tion to the ansatz ([5]) it is worthwhile to mention that 
the equation of state of this type which implies that the 
matter distribution under consideration is in tension is 



available in literature and hence the matter is known as a 
'false vacuum' or 'degenerate vacuum' or 'p-vacuum' [T^ 

Now, from the metric ([1]) we get A' — 2Ar, v' = 2Bt 
and = e""^'' and substituting these values in equa- 
tions © - ([7]), together with our assumptions as given in 
equations (|5]) and ([S]), we get 
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84p + p) = 2e'^^\A + B). 
Subtracting equation (ITUl) from (ITTt . we get 

Sttp = e"^''' {2A + 2B) - 
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The equation ([T2|) . alongwith ([T]), then provides 
STTprfei = e-^^" [B\^ ~ ABr^ - 3A] 
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Thus the effective energy density {p)eff, effective ra- 
dial pressure {pr)eff and the effective tangential pressure 
{ptjeff are obtained as 



(P)e// = ^ 
(^)e// = ^ 
iPt)eff = ^ 
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Using equations ((T6| - ((T8)) the equation of state (EOS) 
corresponding to radial and transverse directions may be 
written as 



UJrir) 







[e-Ar^ {2A-^) 


+ 7? 



UJt{r) 



-Ar^ 



(B^r'^ + 2B- ABr^ - A) 



(17) 



(18) 



It is interesting to note that the effective energy den- 
sity (peff), effective radial pressure {pr eff) and effective 
tangential pressure {pt eff) are independent of a. 

We also note that 



dpeff 
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and 



dp, 



r eff 



dr 



< 0. 



We impose the following conditions for our anisotropic 
fluid configuration to be physically acceptable: 
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• The density is positive definite and its gradient is 
negative everywhere within the fluid distribution. 

• The radial and tangential pressures are positive def- 
inite and the radial pressure gradient is negative 
definite. 

The above results and Figs. [1][2] are in agreement with 
these conditions. 

Note that, at r = 0, our model provides 
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0, 
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and 



d^Pr eff 

dr^ 



<0, 



which indicate maximality of central density and cen- 
tral pressure. Interestingly, similar to an ordinary matter 
distribution, the bound on the effective EOS in this con- 
struction is given by < cji(r) < 1, (see Fig. [3]) despite 
the fact that star is constituted by the combination of 
ordinary matter and dark energy. 

The parameter A = j. (pt ef f — Pr ef f) representing 
the 'force' due to the local anisotropy is obtained as 



A 



1 

47rr 



B^r^ - ABr^ - 3A) 



1 

8^ 



(19) 



This 'force' will be directed outward when Pt > i.e. 
A > 0, and inward if Pt < Pr i.e. A < 0. As it is 
apparent from the Fig. 2] of our model with a repulsive 
'anisotropic' force (A > 0) allows the construction of 
more massive distributions. 



By solving Eqs. (|2Tl)-(|23l), we have 
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We also impose the boundary conditions that at the 
boundary {pr)eff{T = R) — and p ef f{r = 0) = b {= a 
constant), where b is the central density. Thus, 
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Combining, equations p4| and p7)) . we get 
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Note that the values of B obtained from equations ((25)) 
and (|28p are identical. 

At this juncture, to get an insight of our model, let us 
first evaluate some reasonable set for values of A, B, and 
C. According to Buchdahl |18|, the maximum allowable 
compactness (mass-radius ratio) for a fluid sphere is given 
by < |, Accordingly, let us assume that we have a 
dark energy star whose mass and radius are such that 
^ = 0.3999052. Due to highly compact nature of the 
star, we set the radius of the star at i? = 8 km. With 
these specifications, we obtain the values of the constants 
A, B and 6 as ^ = .025, B = .030883, b = .002984. Later, 
we have shown that these values of A and B are justified 
since the energy conditions imply 2A > B > (see Sec. 
V). 



III. BOUNDARY CONDITIONS 



IV. TOV EQUATION 



We match the interior metric to the Schwarzschild ex- 
terior 



ds^ 
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at the boundary r = R. Assuming continuity of the 
metric functions gtt, Qrr and at the boundary surface 
5, we get 
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For an anisotropic fluid distribution, the generalized 
TOV equation is given by 

+ ^' {Peff + Pr eff) + ^ {Pr eff " Pt eff) = 0. 

_ (30) 
Following Ponce de Leon [1^, we write the above TOV 
equation as 

^lajPeff +Pr eff) _ dpr eff 

2 

+ -iPteff-Pr eff)^0, (31) 

where Mq — MQ{r) is the effective gravitational mass 
inside a sphere of radius r and is given by 



1 2 ^ 



(32) 





FIG. 2: Effective radial pressure and transverse pressures are 
plotted against r. 



FIG. 4: The variation of the force, A = | (pt e// — Pr eff) 
due to the local anisotropy with respect to r. 



where, 



which can easily be derived from the Tolman-Whittaker 
formula and the Einstein's field equations. Obviously, the 
modified TOV equation describes the equilibrium condi- 
tion for the dark star subject to gravitational and hydro- 
static plus another force due to the anisotropic nature of 
the stellar object. Using equations (ITU) - ([TE|. the above 
equation can be written as 



F, + Fh + Fa = 0, 



(33) 



= 

Fh = 

Fa - 



-Br {peff+Pr eff); 
dPr eff 

dr ' 



(34) 
(35) 

(36) 



The profiles of Fg, Fh and Fa for our chosen source are 
shown in Fig. [S] The figure indicates that the static equi- 
librium is attainable due to pressure anisotropy, gravita- 
tional and hydrostatic forces. 
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FIG. 5: Three different forces acting on fluid elements in static FIG. 6: The variation of left hand side of the expressions of 
equilibrium is shown against r. energy conditions are shown against r. 



V. ENERGY CONDITIONS 

In this section, we verify whether our particular choices 
of the values of mass and radius leading to solutions for 
the unknownn parameters, satisfy the following condi- 
tions through out the configuration: 

Peff > 0, 
Peff +Pr eff > 0, 
Peff+Pt eff > 0, 
Peff + Pr eff + 2pt eff > 0, 



Employing the energy conditions at the centre {r — 0), 
we may get a bound on the constants A and B as follows: 

(i) NEC: Poeff+Poeff>O^A + B>0, 

(ii) WEC: Poeff+PQeff>O^A + B>0, 

POeff>0^ A>0, 

(in) SEC: po eff+Poeff>O^A + B> 0, 

3po eff + PQ eff >0 ^ B >0, 

(iv) DEC: Poeff> \Po eff\^2A>B. 



Peff > \Pr e//|, 
Peff > \Pt eff\. 

Note that all the energy conditions namely, the null en- 
ergy condition (NEC), weak energy condition (WEC), 
strong energy condition (SEC) and dominant energy con- 
dition (DEC), for our particular choices of the values of 
mass and radius, are satisfied as shown in Fig.|6l It is in- 
teresting to note here that the model satisfies the strong 
energy condition, which implies that the space-time does 
contain a black hole region. 

The anisotropy, as expected, vanishes at the centre i.e., 
Pt eff = Pr eff = Po ef f = at r=0. The effective 

energy density and the two pressures are also well be- 
haved in the interior of the stellar configuration. 



VI. STABILITY 

For a physically acceptable model, one expects that the 
velocity of sound should be within the range < Vs = 
(^) < 1 [13, mi- In our anisotropic model, we define 
sound speeds as 

2 ^ dpr eff 

dpeff 

1 , AAre-^^\A + B) 

{2A- ^)2Are-A-' + ^ {I - e-A-y ^ ' 
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2 _ dpt eff 
dpeff 

e"^'^' [2Ar {B^r^ + 2B ~ ABr^ - ^) + 2Br{A - B)] 

" {2A - ^) 2Are-^-' + ^ (l " e-^''') ' 

(38) 

We plot the radial and transverse sound speeds in Fig. [7] 
and observe that these parameters satisfy the inequalities 
< ^ 1 ^nd < < 1 everywhere within the stellar 
object. 

Equations dHT]) and lead to 

2 2 

e"'^''' [2A^Br^ + GAV + 2Sr2 - 2ABr - 2AB^r^] 

{2A - ^) 2Are-A-' + ;! (l " ^-^H ' 

(39) 

From equation (|39)) . we note that v'^f — v-^j. < I . Since, 
< ^ 1 and < ii^t < 1, therefore, | — |< 1. 
In Fig. m we have plotted | — ^sr I- 

Now, to examine the stability of local anisotropic mat- 
ter distribution, we use Herrera's [20| cracking (or over- 
turning) concept which states that the region for which 
radial speed of sound is greater than the transverse speed 
of sound is a potentially stable region. Thus, if the dif- 
ference of the two sound speeds -y^i — retains the same 
sign everywhere within a matter distribution, no crack- 
ing will occur. In our case. Fig. |9] indicates that there is 
no change of sign for the term — v'^^, within the specific 
configuration since the difference is negative throughout 
the distribution. Therefore, we conclude that our dark 
energy star model is stable. 

VII. MASS-RADIUS RELATION 

In this section, we study the maximum allowable 
mass-radius ratio in our model. For a static spherically 
symmetric perfect fluid star, Buchdahl [31 showed that 
the maximally allowable mass-radius ratio is given by 
< I (for a more generalized expression for the same 
see Ref. [23). In our model, the effective gravitational 
mass in terms of the effective energy density pe/ / can be 
expressed as 

M,ff ^AnJ (p + pde) r^dr = - e"^"" J • 

(40) 

In Fig. [TUl we plot this mass-radius relation. We have 
also plotted against R (see Fig. [TT|) which shows 

that the ratio —j^ is an increasing function of the radial 
parameter. We note that a constraint on the maximum 
allowed mass-radius ratio in our case is similar to the 
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FIG. 7: The variation of radial sound speed v^,. and tangential 
sound speed v^t are shown against r. 
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FIG. 8: The variation of j v^t — f ?r I is shown against r. 

isotropic fluid sphere, i.e., ^ < | as obtained earlier. 
The compactness of the star is given by 

u=^^ = l{l-e---'). (41) 

The surface redshift (Zg) corresponding to the above 
compactness (m) is obtained as 

= (1 - 2u)-^ - 1, (42) 

where 

Zs = e4«^ (43) 
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FIG. 10: The variation of M^ffecUve is shown against R. 

Thus, the maximum surface redshift for an anisotropic 
star of radius 8 km turns out to be Zg — 2.225541. 



VIII. JUNCTION CONDITION 

One of the issues in connection with a static anisotropic 
matter distribution is that, though the radial pressure 
at the boundary of the star must vanish, the tangential 
pressure is not necessarily zero at the boundary. This 
forces us to examine the junction conditions of a static 
anisotropic star in closer details. We propose here a shell 
type envelope at the boundary surface so as to address 



FIG. 12: The variation of redshift function Zs is shown 
against R. 



this issue. Note that the fundamental junction condi- 
tion for a static star is that there has to be a smooth 
matching between the interior solution and Schwarzschild 
exterior at the boundary. Now, though the metric co- 
efficients must be continuous at the junction surface S 
where r = R, their derivatives may not be continuous 
at the junction. In other words, the affine connections 
may be discontinuous at the boundary surface. This can 
be taken care of if we consider the second fundamental 
forms of the boundary shell. The second fundamental 
forms associated with the two sides of the shell [23l - [26j 
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are given by 



a/3' 



(44) 



Is 



where are the unit normals to S and can be writtten 
as 



df df 



df 

dX'^' 



(45) 



with n^'n^ = 1. In Eq. ps)) . i^' are the intrinsic coor- 
dinates on the sheh with f — being the parametric 
equation of the shell S and — and + correspond to in- 
terior and exterior (Schwarzschild) metrices. Note that 
radial pressure on the shell is zero. By using Lanczos 
equations (23l - [26| , the surface energy term E and surface 
tangential pressures pe = pcf, = pt may be obtained as 
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Since the metric functions are continuous on S, we have 

E = 0, (48) 

and 



1 



Pt 



1 _ 2M 
^ R 



(49) 



Therefore, we are now in a position to match our interior 
solution to the Schwarzschild exterior in the presence of 
a thin shell. 

From Fig. [T31 we note that the effective transverse 
pressure at the boundary (i? = 8 km) is positive though 
the transverse pressure in the shell is negative. This 
clearly indicates that the static equilibrium may be at- 
tained due to positive pt eff and negative pt shell- In 
this figure, we have plotted pt eff and pt shell together 
by assuming that the width of the shell is in between 
8—15 km. We observe that at r = 15 km, both pt eff 
and Pt shell vanish simultaneously. Thus the thickness of 
the shell in this case is 7 km. Physically, this suggests 
that anisotropic matter is confined within 8 km from the 
centre of the star and the outer region contains a thick 
shell extending up to 7 km. The thick shell is charac- 
terized by zero energy density and non zero transverse 
pressure though the shell does not exert any radial pres- 
sure. 



IX. CONCLUSION 

Dark energy stellar models have found astrophysical 
relevance for various reasons, one particular reason be- 
ing it's importance as an alternative candidate to a black 
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FIG. 13; Plot for pt eff and p 



is shown against r. 



hole. The model developed here satisfies all the physi- 
cal requirements and is horizon free and, therefore, can 
potentially describe a compact object which is neither a 
neutron stars nor a quark star. 

Similar to many model proposed earlier, our model also 
requires an envelope just outside the baryonic matter for 
smooth matching with the Schwarzschild exterior space- 
time. For matching we have firstly, assumed continuity 
of the metric functions gtt , grr and at the boundary 
surface S*, and secondly imposed the boundary conditions 
that at the boundary {pr)eff{'r — R) = and p eff{r — 
0) — b {= a constant), where b is the central density. 
Thus we get a set of expressions for A, B and C which, 
due to the maximum allowable compactness for a fiuid 
sphere 'l8], eventually can be worked out as A = .025, 
B = .030883, b = .002984 for the assumed mass-radius 
ratio as ^ = 0.3999052. Later on, in Sec. V, it has been 
shown that these values of A and B are justified since 
the energy conditions imply 2 A > B > 0. 

Regarding stability of local anisotropic matter distri- 
bution, we use cracking concept of Herrera [20| which 
states that the region for which radial speed of sound is 
greater than the transverse speed of sound is a poten- 
tially stable region. It is observed from Fig. 12] that there 
is no change of sign for the term v'^^ ~ '^sr within the 
specific configuration and hence advocating in favour of 
stability of our dark energy star model. 

Let us now concentrate on some of the other works on 
KB analysis, especially the works by Varela et al. [27} and 
Farook et al. In both the works static, spherically 

symmetric, Einstein-Maxwell spacetime have been con- 
sidered with a fluid source of anisotropic stresses whereas 
the present investigation is neutral one with anisotropic 
fluid source. However, a common feature of all these KB- 
models is singularity-free, stable conflgurations. Valera 
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et al. 27[ in their work have found out a link of their con- 
struction with a charged strange quark star as well as 
models of dark matter including massive charged parti- 
cles. Farook et al. 28], by using a Chaplygin-type EOS, 
predicted the possible existence of a Chaplygin charged 
dark energy star or a strange quark star of radius about 
8 km. 

Therefore, it is interesting to note that present work 
deals with a singularity-free spherically symmetric body 
of radius r = 15 km such that both pt eff and pt shell 
vanish simultaneously. The model physically contains 
anisotropic matter which is confined within 8 km from 
the centre of the star and the outer region contains a 
thick shell extending up to 7 km. Here the thick shell is 
characterized by zero energy density and non-zero trans- 
verse pressure though the shell does not exert any radial 
pressure. Note that in a recently proposed toy model 
[29| , a relativistic stellar configuration has been devel- 
oped where the core of the star is characterized by a 
wormhole like solution for some kind of exotic matter vi- 
olating the weak/null energy condition and is surrounded 



by some ordinary matter satisfying a polytropic EOS. 
This kind of theoretical modelling would get observa- 
tional support in the future. 

We hope our model inspires observational workers to 
search this type of stars. That is, we mean, the stars 
containing anisotropic matter which is confined within 
certain radius from the centre of the star and the outer 
region contains a thick shell extending up to several kilo- 
meter. 
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